Quantum fluctuations and coherence in high-precision single-electron capture 
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The phase of a single quantum state is undefined unless the history of its creation provides a 
reference point. Thus quantum interference may seem hardly relevant for the design of determinis- 
tic single-electron sources which strive to isolate individual charge carriers quickly and completely. 
We provide a counterexample by analyzing the non-adiabatic separation of a localized quantum 
state from a Fermi sea due to a closing tunnel barrier. We identify the relevant energy scales and 
suggest ways to separate the contributions of quantum non-adiabatic excitation and backtunneling 
to the rare non-capture events. In the optimal regime of balanced decay and non-adiabaticity, our 
simple electron trap turns into a single-lead Landau-Zener-backtunneling interferometer, revealing 
the dynamical phase accumulated between the particle capture and leakage. The predicted "quan- 
tum beats in backtunneling" may turn the error of a single-electron source into a valuable signal 
revealing essentially non-adiabatic energy scales of a dynamic quantum dot. 

PACS numbers: 73.63.Kv, 73.23.Hk, 73. 21. La 



Successful demonstration of electron-on-demand 
sources based on electrostatic modulation of nano- 
electronic circuit elements such as dynamic quantum 
dots [TH3] or mesoscopic capacitors [H [5] has offered a 
prospect of building an electronic analog of few-photon 
quantum optics [6] that exploits the particle-wave 
duality and entanglement of individual elementary 
excitations in a Fermi sea [TrfTU]. This ambitious goal is 
complemented by a long-standing challenge in quantum 
metrology pd to untie the definition of ampere from 
the mechanical units of SI [12] and implement a current 
standard based on direct counting of discrete charge 
carriers. So far the overlap between these research direc- 
tions [331 [14] has been rather limited arguably because 
metrological applications strive to maximize the particle 
nature of on-demand excitations. Optimizing the 
trade-off between speed and accuracy of single-electron 
isolation [2] does require consideration of quantum error 
mechanisms such as non-adiabatic excitation j!5H17j or 
backtunneling [18-21 ]. However, these effects have been 
hard to differentiate experimentally owing to complexity 
of non-equilibrium many-particle quantum dynamics 
[22) and experimental challenges in exercising high-speed 
control of the electrostatic landscape. The quantum 
phase of the captured particle has been considered thus 
far as inconsequential for accuracy and inaccessible 
for measurement unless the particle is ejected into a 
separate interferometer [TTJ] . 

In this work, we propose a new type of interferometry 
to measure and thus control the non-equilibrium energy 
scales governing the decoupling of a dynamic quantum 
dot from a Fermi sea. Remarkably, our approach requires 
neither multiple spatial paths [TU] nor noise measure- 
ments [H [22 [23], relying instead on quantum beats in 
spontaneous emission of electrons back to the source lead. 
Using a generic (35[ [53] effective single-particle model 
we predict an interference pattern in the charge capture 
probability that reflects the dynamical quantum phase 
accumulated between the isolation of the localized quan- 
tum state and the onset of backtunelling. 
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FIG. 1. (color online) (a) Schematic plot of the rising 
level energy e(t) and exponentially decreasing broadening 
T(t). Pictograms illustrate the two paths for interference: 
dark (blue) arrows mark non-adiabatic excitation followed 
by phase accumulation above the Fermi edge /i, light (yel- 
low) arrows show adiabatic "elevator ride" and backtunnel- 
ing. (b) Instantaneous current flowing into the quantum dot, 
h(t). Thick (blue) line - exact Eq. p^l, thin dashed line - 
quantum-broadened Markov approximation, Eq. Q. Model 
parameters: e c /A pt b = — 6, r c /A pt b = 1.2, T = 0. 



Our innovations can be explained within a simple two- 
path picture of Mach-Zchnder interferometry in time do- 
main, see Fig. [TJa). The first branch is quantum ex- 
citation above the Fermi edge in the source lead due 
to finite time-scale r for the pinch-off of tunneling [TB] 
[see dark (blue) arrows in Fig. [it a)]. The correspond- 
ing energy spread and the (small) path-splitting ampli- 
tude can be estimated by the energy-time uncertainty, 
T c = %/{-kt), and the Landau-Zener theory [35], respec- 
tively. The second branch is adiabatic lifting of the oc- 
cupied energy level followed by splitting of a small am- 
plitude back into the lead once the level emerges above 
the Fermi sea [see light (yellow) arrows in Fig. [It a)]. 
This elevation-backtunneling branch is sensitive to elec- 
trostatic cross-coupling between the barrier and the bot- 
tom of the confining potential well; we quantify this 
crosstalk by plunger-to-barrier ratio A pt b = st (shift 
of the localized energy level e during the characteristic 
decoupling time). The output ports of the interferom- 
eter can be read either by detecting excitations created 
in the Fermi sea, or, more conveniently, by measuring 
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the charge capture probability. The latter is accessi- 
ble in experiment by repeated ejection of the captured 
electrons into a collector lead and measuring the result- 
ing dc current [Tl |21 [Til ITM2Tj . The beam-splitters are 
tuned by A pt b/r c to maximize contrast, the phase mea- 
surement reveals the "elevator speed" i oc A pt br c , and 
the temperature smearing allows absolute calibration of 
the energy scales. Our scheme is conceptually related to 
Landau-Zener-Stiickelberg interferometry [2"6Tf2"5] which 
measures the relative dynamical phase of discrete states 
via creation of superposition in sequential non-adiabatic 
level crossings. In contrast, we propose to access the 
phase of a single localized state measured against a ref- 
erence point in the continuum (defined by a sufficiently 
fast decoupling and a sharp Fermi edge). 

The proposed measurement addresses a persistent 
challenge in robust utilization of electrostatically defined 
quantum dot devices: control of rate (r), type (barrier 
versus plunger) and magnitude (A pt b) at which exter- 
nal voltage pulses are converted into the time-dependent 
potential guiding individual transport electrons on the 
chip. Although parametric time-dependency of the elec- 
tronic matrix elements is the standard input for theory 
[2H HI] , in practice mesoscopic fluctuations (due to dis- 
tribution and charge switching dynamics of impurities, 
finite fabrication precision etc.) and challenges of signal 
propagation at high frequencies (GHz range) often re- 
quire measuring the characteristic quantities on sample- 
to-sample basis. For near-equilibrium, bias spectroscopy 
[50] provides a versatile tool, but for the large-amplitude, 
high-frequency modulation the options are limited [T]. 
As a foreseeable direct application of our results we ex- 
pect a reliable measurement of r c /A ptb to help settling 
the debate on the fundamental factors limiting the preci- 
sion of the state-of-the-art single-electron-based current 
sources [2]. More generally, we hope our analysis will fa- 
cilitate the crossover of ideas between fundamental and 
applied directions of single electronics. 

Model and formalism. The model is described by 
an effective single-particle Hamiltonian % = e(t)d'd + 

J2ki e k c k c k + ^M[ c !^ + ^ c fe]}' where S creates a lo- 
calized non-degenerate electronic state in the dot and 
c k — a quasi-continuous state in the lead. The lead 
is connected to a thermal reservoir with chemical po- 
tential fi — and temperature T. Employing a time- 
dependent tunneling Hamiltonian relies on timescale sep- 
aration [21] : fast screening in the leads determines the in- 
stantaneous values of the slowly varying parameters for 
the underscreened region (the quantum dot). We choose 
/c-independent real V and the wide-band limit so that 
r(t) = 2irpV 2 (t) and e{t) are the fully dressed elastic 
width and the on-site energy respectively (p is the den- 
sity of states in the lead) . 

The quantum kinetic equation for the average occupa- 
tion within the model, n(t) = (cv (t)d(t)) , is given by the 



non-equilibrium Green functions theory [24] . 
fin{t) = -T{t)n(t)- 

Im f JT(t)T(t')G(t, i')e ie( *-*' )/ft dt'de , (1) 

where /(e) is the Fermi distribution, and G(t, t') is the 
retarded Green function of the level, 



)/2]/fi 



G{t,t') = _ie(t-i / )e~ iJ '*' <ttlKtl) ~ <r{ * 1 
Integrating Eq. (I} gives [cf. Eq. (44) of Ref. 



(2) 



n(t) = 



dt' v / T(f)G(t,t r )e^ et ' /h de. (3) 



Equation ^ is the sum of probabilities to be scattered 
into the localized state from an occupied state in the 
continuum; one can show [5T][52] (see Appendix [Aj that 
it agrees with the Floquet formalism [29] [33] which is 
often used in the scattering form to study single-charge 
emitters [25] [5H [55] . In this work we consider electron 
trapping achieved by reducing T(t) to zero as t — > oo and 
compute the capture probability rif = lim t _ i . 00 n(t). 

The history of parametric time-dependence that af- 
fects nit) is limited by a finite memory time r mem = 
ftmm{(fcT) _1 ,r _1 (t)} due to (a) lead-induced dephas- 
ing of the discrete state, and (b) thermal smearing in the 
reservoir. One can show (see Appendix [b]) that for suffi- 
ciently slow processes, when |e|T^ em -C fi (no phase rota- 
tion) and r|-r mom <C T (well-defined T), the exact kinetic 
equation (TJ) can be replaced by a (quantum broadened) 
Markov approximation, 



hn{t) = -r(t) {n(t) - n cq [e(t),T(t)]} , 



(4) 



with parametrically defined standard equilibrium occu- 
pation, n cq {e,T) = (27T)" 1 / dw Tf(u)/[(u - e) 2 + T 2 /4]. 
Although Eq. p} still permits strongly non-adiabatic sce- 
narios [36] . for the decoupling problem at hand r mcm di- 
verges at zero temperature as t — > oo, thus we must use 
Eq. |3j. 

We define a crossover moment t c as the earliest time 
from which the quantum phase of the level can be pre- 
served, J. °°T(t)dt = h, and explore exponential time- 
dependence [16] of T(t) around t c , 



r(t) = 7rT c 



a -(t-te)A 



(5) 



accompanied by a linear shift of energy [TH] [2H1 [35] [SS] , 



e(t) =e c + A pth (t~t c )/i 



(6) 



The shape of e(t) and T(t) for t <C t c as well as the initial 
conditions for Eq. ([T]) are irrelevant for nf if the ansatz 
^ holds from a few t before t c . Whether the last particle 
exchange between the dot and the source (most likely 
to occur around t ps t c ) results in a captured electron 
(rif ~ 1) or a hole (n* w 0) depends on the position of 
e c with respect to \i (below or above, respectively). 
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Experimental realization. A prototypical realization 
for the model is a small near-empty electrostatically de- 
fined quantum dot with large level spacing and charg- 
ing energy, as in, e.g. [31 [T7]. A linear ramp of volt- 
age Vi(i) on a gate that defines the tunnel barrier be- 
tween the source and the quantum dot creates the time- 
dependencies ([5]) and (|6| while a static voltage V 2 on 
another gate can be used to tune the decoupling energy, 
e c = — aV^Aptb + const. r c ~ t _1 is inversely propor- 
tional to the rise-time of the Vi(t) pulse. 

Once r(i) becomes negligible (t — > 00 in our model), 
further modulation of the confining potential can ensure 
complete ejection of the captured electrons into the drain 
lead [21 El] . Repeating the whole cycle with a frequency 
/ <C t^ 1 and measuring the dc component of the source- 
drain pumping current /(V2) can provide accurate data 
on rif(e c ) = I/(ef) (here e is the electron charge). Note 
that a is the fitting parameter El EQl [37] of the decay 
cascade model [21], see Eq. ( [8c| below. 

Qualitative picture. The essential features of the 
model can be seen in the time-dependent ensemble- 
average current h(t) for Aptb ~ r c an d e c < — A pt b at 
low temperature, see Fig. [lib). Before the formation of 
the localized level, at times t t c — rln |e c |/r c , the av- 
erage occupation number remains constant and equal to 
the average density of electrons per quantum state in the 
lead (1/2 in our dispersionless model). A large peak in 
the current near t c marks adiabatic filling of the rapidly 
narrowing level. Much smaller opposite sign feature indi- 
cates the onset of backtunneling at t > tf, with t^ defined 
by s(tb) = /i. It is instructive to contrast the exact fi(t) 
with the solution of the Markov equation Q shown by 
the dashed line in Fig.Qb). The level of agreement cor- 
relates with the condition T(t) > T c being fulfilled expo- 
nentially well for t <C t c , breaking down around t c , and 
being strongly violated for t > tf,. 

Probability of charge capture. Our main result is 



n f 



dc 



2tt 2 T, 



(7) 



exp 



27rr, 



A 



ptb 



dx 



n/(e c ) is a step-like function changing from 1 to as e c 
goes from —00 to +00. The limit forms are 

|7(e c ), Ap tb ,rW0 (8a) 

n f = J (2/vr)tan- 1 e^/ r =, fcT,A ptb ^0 (8b) 

fcT,rwO (8c) 



The limit (8a) corresponds to a sudden decoupling from 
equilibrium at s(t c ) = £ c , the corresponding step is sym- 
metric under e c — > — e c and maps out the thermal dis- 
tribution [35]. The limit (8b) reproduces the result of 

has been predicted previ- 



Flensberg et al. [16], a symmetric step of width r c . The 
double-exponential shape 



ously [TM2"T] from a master equation and validated ex- 
perimentally [H [37] with up to 10~ 6 relative accuracy. 




FIG. 2. (color online) Charge capture probability nj as 
a function of the level energy e c at the decoupling moment. 
T = and r c /A pt b is varied as indicated. The inset show the 
same quantity in the logarithmic scale with additional thick 
(red) lines marking the asymptotics Q. 



Our derivation reveals the plunger-to-barrier ratio A pt b 
as the energetic measure of the step width. 

Finite temperature is accounted for by simple thermal 
smearing, n/(e c ) = J rif(e c — €)\t=o(— df/de)de, thus we 
focus mainly on T = 0. Figure [2] shows evolution of the 
lineshape nt (e c ) as r c /A pt b is increased. The asymmetry 
around e c = gets "inverted" with respect to the double- 
exponential ([8c]) at r e = (2. . .3)A p tb before approaching 



the symmetric limit ( 8b ) . Non-perturbative asymptotics 
of Eq. Q for e c S> max[A pt b,r c ] gives an exponential 
dependence with a power-law prefactor, 



n f 



2ec 



A p tb/r c 



-e c /r c 



T -> , (9) 



shown in the log-scale inset of Fig. [2] This is similar 
to the Fermi function tail, ~ e _£c / fc , which dominates 
the small-n/ asymptotics for kT > r c regardless of A pt b- 
Thus we conclude that both quantum and thermal fluctu- 
ations always trump the double-exponentially suppressed 
backtunneling for small n/. This finding may have im- 
plications for the minimal slope position on the plateaux 
between successive current quantization steps in single- 
gate pumps [201 EH- A change in slope of log nf versus e c 
as the barrier-closing timescale r is reduced may signal 
that T c /kT > 1, although this effect on its own would be 
hard to differentiate from local heating. 

Quantum oscillations. The asymptotics of 1 — n/ at 
large negative e c switches from ~ e Ec/A pt b ^ ^ e £ c /r c ^ 
T c = A pt b via a surprising sequence of miniplateaus. The 
latter can be seen as ripples in the derivative dnf/de c , 
as shown in Fig. [3](a)-(b). Although the oscillation am- 
plitude decays exponentially, the non-monotonic behav- 
ior of drif(e c )/de c is manifest in the range of T c = 
(0.8 ... 2.2) A p tb at temperatures up to kT w 0.35A pt b, 
see Fig. §c). 

We interpret the oscillations as interference between 
two excitation paths that promote lead electrons above 
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FIG. 3. (color online) (a) The derivative — Aptbdrif /de c 
(vertical axis, log scale) as a function of the decoupling energy 
e c and the ratio r c /A p tb at T = 0. Three cuts at r c /A p tb = 
0.9, 1.4 and 2.1 are shown on a linear scale in (b). (c) The ratio 
of the first interference maximum in —dnf(e c )/de c to the first 
minimum on a log scale as a function of r c /A pt b as kT/A pt b 
increases through 0, 0.1, 0.2 and 0.3 (from the topmost to the 
lowermost curve). 
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FIG. 4. Main panel: the derivative —dnf/de c multiplied by 
A p tbe" £c/rc for r c /A ptb = 1.4 and T=0. Left inset: the po- 
sitions Em of the minima in the main graph versus their index 
number, m, scaled according to the two-path formula (13 1. 



Right inset: the absolute difference between a maximum and 
the preceding minimum of the oscillations shown in the main 
panel on a log scale versus a scaled — e m . The straight line of 
slope of ±1 in both insets is plotted according to Eq. (131. 



the Fermi energy. One path is capture-elevation- 
backtunneling similar to elevator resonance activation 
proposed by Azbel' [33]. As seen from the scattering 
interpretation of Eq. ([7]), the electrons most likely to be 
captured have incoming energies e ~ e c , thus we estimate 
the amplitude to be raised by an "elevator ride" from e c 
to an energy e e > n by a three-amplitude product, 

^eicv(e e ) oc V(t c )e- 1 ^ e{t)dt/n V{t e ) , (10) 

where the exit time t e > tb is determined from e(t e ) = e e . 

The other path from e ss e c to e = e e is Landau- 
Zener-like excitation due to time-dependence of T(t). We 
estimate the corresponding amplitude by following the 
Landau solution of a two-level problem [25 with ma- 
trix elements H\\ — e e , H 22 = s c , and H 12 = H 2 i = 
r(t) [IS]- The adiabatic eigenvalues E\(t), E 2 (t) have 
a gap AE(t) = ^4T 2 (t) + (e e -e c ) 2 > on the real axis 
but become degenerate if analytically continued into the 
complex t-plane. The branching point of E(t) with the 
smallest positive imaginary part, to = t c + r ln[2r c /(e e — 
s c )] + iirr/2, determines the transition probability [25] 
exp[-2Im J^ ctg AE(t)dt/h] = exp[-7rr(e e -e c )/K\. Ne- 
glecting pre-exponential and logarithmic terms, the am- 
plitude for excitation at t c and evolution up to t e is 

V L z(e e ) cx e -(^-^)/(2r c )-^o e -ie e (t e -t c )/n ^ (n) 

where ipo is the phase of the Landau-Zener transition 
(Stokes phase) known to depend weakly on energy [2"5] , 
We estimate the total non-capture probability as 

/>oo 

1-n/oc / de e \ip LZ {e e ) + ip c i cv {e e )\ 2 - (12) 

The competition of |t/>lz| 2 and IV'oicvl 2 at large nega- 
tive e c agrees with the asymptotic envelope of nt while 



R' e V'Lz(0)^/'eiev('-') gi yes t ne oscillating part of dnf/de c , 

e-^ +r --'" 2 cc S (, + ^L-). (13) 

The argument of the cosine is the dynamical phase ac- 
cumulated between t c and tb, see the shaded triangle in 
Fig. [T] Despite the crudeness of approximations leading 
to Eq. p"3] ) , the result agrees well with the exact solution 
Eq. Q both in amplitude and phase, as shown in Fig. [4] 

Feasibility. The oscillation effect does not rely on 
any fmetuning (apart from balancing the interferometer, 
Aptb ~ F c ) thus it should be robust against deviations 
from the assumed time-dependencies. The most impor- 
tant foreseeable limitation is the overlap of additional 
interference modes which must be separated by a suffi- 
ciently large on-the-dot level spacing Ae from t t c on- 
ward; our single-mode formula Q is limited to — e c < Ae 
requiring Ae > 8A pt h to resolve the first interference 
minimum (see AppendixjC]). Using a recent report [17] on 
out-of-equilibrium excited states in a rapidly decoupling 
dynamic quantum dot we read a = 0.28 mV _1 from 100 
MHz data in Fig. 1 of Ref. ITTland estimate the gap to the 
first exited state observed at 1 GHz to be Ae/A pt b ~ 4.5 
which is the same order of magnitude as our require- 
ments. Measurements of a similar device [5] have just 
reached the required accuracy threshold for n/(e c ) sug- 
gesting the proposed dynamical phase interferometry is 
feasible with current technology. 

Conclusions. We have considered the quantum dy- 
namics of isolating a single particle from a Fermi sea by 
a closing tunnel barrier and proposed ways to measure 
the relevant non-equilibrium energy scales. Time-domain 
interferometry revealing the dynamical phase of the con- 
fined particle offers a unexpected "built-in" instrument 
to quantify and separate the quantum effects hampering 
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deterministic electron-on-demand operation which is a 
cornerstone for future on-the-chip electron quantum op- 
tics and a quantum realization of the ampere. 
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Appendix A: Relation to Floquet formalism 

Here we demonstrate the connection of Eqs. (1) and (3) 
in our paper to the instantaneous current formulae de- 
rived in Ref. [23] for a periodically modulated mesoscopic 
system. Using Eqs. (9) and (10) of Ref. [531 the average 
incoming current can be expressed as 



(/(*)> 



e 
2irh 



de 



U m (e — nihil) 



1 



(Al) 



where 27r/fi is the period of parameter modulation. The 
Floquet scattering matrix U m (e) is a double- Fourier rep- 
resentation of the time evolution operator for the lead 
electrons [see Eqs. (5)-(6) of Ref. [23], 



U(t,t') 



/rfwe- ia ' (t -*') +imm 'U m (M- (A2) 



On the other hand, scattering amplitude in the lead and 
the retarded Green function on the quantum dot are re- 
lated by the standard equation 

U{t, t') = 8{t - t') - 2iri P V(t)G(t, t')V{t')/h , (A3) 

which follows directly form the equation of motion for 



the Heisenberg operators in the leads. Using Eqs. (A2) 
and ( A3 1 to expresses the time-dependent T-matrix 



T{e,t) = TT x jf(i) 



G(M')V , rMe K(i "'' )/ ' , 'it' 



■lOI 



U m (e — mhfl) , 



(A4) 



allows re-writing Eq. (|A1|) as 

(/(*)) 



e 

2nh 



def(e) [|TM)| 2 + 2ImTM)] 



(Al) 



Since our result (3) for the time-dependent level occupa- 
tion can be written as 



T(t)n(t) 



1 

2^ 



/(e)|TM)| 2 de, 



(3') 



equivalence between the Floquet current (Al) derived in 
Ref. [23] and the quantum kinetic equation (1) from the 
non-equilibrium Green functions theory is now explicit, 
<J(i)> = -en(t). 



Appendix B: Markov approximation 

The role of quantum effects becomes clearer if we con- 
trast Eq. ([!]) with a Markov master equation. Performing 
the frequency integral in Eq. (1) one gets 

fin = -T(t)n(t) + 2 Re / dt'y/T(t)T(t')f*(t - t') 



xexpj-i J dt[e(t) -iT{t)/2]/K\ , (1') 

where f(t) = S(t)/2 + ikT/[2hsmh(irkTt/h)}. The main 
contribution to the memory integral in Eq. ( [T 7 ] ) is from 
an interval of length r mcm = 7Lmin{(fcT) -1 , r _1 (i)}. If 
the time-dependence of e(t') and T(t') can be ignored 
around t « if over a time span longer than T mem then 
the exponential in Eq. ( [T 7 ] ) can be approximated by 
e -i[e(t)(t-t )+T(t)\t-t |/2]/B anc j memor y integral eval- 
uated as if the time-dependences were frozen. This pro- 
cedure results in the quantum broadened Markov ap- 
proximation, see Eq. Q of the main text. The neces- 
sary conditions are |e|T^ em <C h (no phase rotation) and 
|r|r me m r (well-defined T), or, for our special model, 

r c , (reAptb) 1 / 2 « max[r(t), kT] . (Al) 

Thus the Markov approximation Q is adequate when 
the quantum broadening due to time-dependence of rates 
(h/r) and energy (-y/?i|e|) is smaller than the quasistatic 
level width (thermal, kT, or tunneling, T). 



Appendix C: Effect of excited states 

Proper investigation of the interplay between the dy- 
namical interference effect discussed in the paper and the 
exited states in the quantum dot is a open research topic 
beyond the scope of the present communication. Never- 
theless, we can make a simple estimate of the effect of 
finite Ae in the regime when backtunneling dominates 
(r c -C Aptb) using the standard Markovian master equa- 
tions [41] . 

We consider two levels £\{t) and £2(^)1 separated by a 
constant gap Ae, 

e 1 (t)=e c + A pth (t-t c )/T, (Ala) 
62 (*) =£c + Ae + Aptb(t-*c)/r. (Alb) 

We assume equal tunnel broadenings for the levels, each 
equal to T(t) as specified in Eq. (4). In the limit of large 
charging energy, the kinetic equations |41) for the indi- 
vidual level occupation probabilities are (i = 1,2): 

&P i (t)=r i (t)P i (f)[l-/(e i (t))] 

+ r i (t)[l-n(t)]/( £i (*)) (A2) 

with n(t) = Pi(t) + P2(t). An initial condition at a suf- 
ficiently early time (we take t = — 20r) can be quanti- 
fied by an effective excitation temperature T* (in energy 
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T"= 1 

ft(to) = 3.4-l(H 




de c 



1.5 



5 Ae = 8 10 



0.5 



0.0 




FIG. Al. Charge capture probability computed within a two-level model defined by Eqs. |5|, \A1\ , and |A2l. All energies are 
in units of A pt b and A pt b 3> T — > 0. 



units), P 2 (t )/Pi(t ) = e- Ae/T *, Pi(t ) + P 2 (t ) = 1. 
No further mechanisms of excitation or relaxation be- 
sides the charge tunneling [Eq. ( |A2[ )] are included in the 
model. 



Solving (A2) numerically for Ae > A pt b results in 
the same double-exponential shape for n/(e c ) as for the 
single-level model, Eq. (10c), except for an additional dip 
at — e c w Ae. The depth of the dip is proportional to the 



initial fraction of the exited state, ^2(^0)7 see Fig. Al '&). 
For this reason we have interpreted the distance between 
the ground state and the first excited state features in 



the charge capture probability measured in Fig. 1(g) of 
Ref. [T7| as a direct indicator of Ae and thus obtained an 
estimate Ae/A pt b ~ 4.5 quoted in the main text. 

We have verified that, on level of master equa- 



tions (A2), the exponentially boosted derivative signal 
e~ Be ' s> tb dn f I de c (analogous to the one used in Fig.Ulto 
demonstrate interference) is not qualitatively affectea by 
the presence of an excited state if — e c < Ae — 3A pt b, see 

c / Aptb I for which the 
5 [see Fig. m 



Fig. Al'b). Since the smallest 



first interference minimum appear is 
we may expect key results of our single-level modeling to 
hold for Ae > 8A pt b, as stated in the main text. 
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